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Abstract: In this work, we present a realistic analysis of the potential of the present- 
day reactor experiments Double Chooz, Daya Bay and RENO for probing the existence of 
sterile neutrinos. We present exclusion regions for sterile oscillation parameters for each 
of these experiments, using simulations with realistic estimates of systematic errors and 
detector resolutions, and compare the sterile parameter sensitivity regions we obtain with 
the existing bounds from other reactor experiments. We find that these experimental set- 
ups give significant bounds on the parameter O ee especially in the low sterile oscillation 
region 0.01 < Am 2 41 < 0.05 eV 2 . These bounds can add to our understanding of the sterile 
neutrino sector since there is still a tension in the allowed regions from different experiments 
for sterile parameters. 
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1 Introduction 

For a long time, it was taken as an accepted fact that there are three families of fermions, 
and hence three neutrino flavours. But experiments like LSND [1] and now MiniBooNE 
[2] have indicated the presence of a fourth type of neutrino - the sterile neutrino. Sterile 
neutrinos do not carry Standard Model gauge quantum numbers, so they do not take part 
in Standard Model gauge interactions, but they can mix with the other 3 active neutrinos. 
The LSND result provided evidence for — > v e oscillations from /i + decay-at-rest (DAR), 
and the same oscillation channel was studied in the MiniBooNE experiment. Similar results 
are also indicated by i> e and v e disappearance channels, by the reactor anomaly [3] and 
Gallium anomaly [4, 5], and from cosmology [6-10]. These results suggest that the deficit 
of observed neutrino fluxes in the respective channels {y e or i7 e ) may be an indication of 
the existence of a fourth type of neutrino. We know that the peak energy of i7 e produced 
during beta decay in reactors is ~ 3.5 MeV. If the typical baseline of a neutrino oscillation 
experiment is ~1 m, then the order of Am 2 which could be detected is ~ ^MeS/ ~ 1 e ^ 2 - 
This is the relevant scale for the oscillation of active neutinos to sterile neutrinos. 

The latest global fits of sterile neutrino parameters were presented in [11, 12]. In [11], a 
fit of all present and past experiments giving bounds on the sterile parameters is performed 
in a (3+1) scenario (i.e., 3 active and 1 sterile neutrino flavors). The short baseline neutrino 
oscillation results from KARMEN [13, 14], LSND [15], v e + 12 C -+ 12 N g . s , + e~, as suggested 
in [16] have also been included. Their best fit values at 95 % CL, including MiniBooNE 
u e and u e data, are sin 2 20 ee = 0.1 and Am| 1 =0.9 eV 2 . Earlier, such global analyses on 
sterile neutrino parameters in (3+1) or (3+2) scenarios have been presented in [17-23]. 
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A significant amount of work is already available in the literature regarding the search 
for sterile neutrinos in reactor and atmospheric neutrino oscillation experiments. Recently 
the possibility of using atmospheric neutrinos as a probe of eV 2 -scale active-sterile oscil- 
lations was explored in [24], where bounds on sin 2 20^ and Am^ were presented. The 
implication of sterile neutrinos on measurements of #13 in a (3 + 2) scenario in the Double 
Chooz [25] reactor experiment was studied in [26]. The impact of light sterile neutrinos 
on measurements in Double Chooz and Daya Bay [27] was studied in [28] in a (3+1) 
scenario. A study of the effect of sterile neutrinos on #23 and #13 measurements in MINOS, 
T2K and Double Chooz was performed in [29]. Similar studies for Daya Bay were carried 
out in [30]. A search for sterile neutrinos using MINOS was done in [31]. A constraint on 
the mixing angle #14 from a combination of Solar and KamLAND data was given in [32, 33]. 
An analysis of the results of Double Chooz, Daya Bay and RENO to simultaneously fit #13 
and the reactor neutrino anomaly was recently performed in [34]. 

In this work, we present exclusion regions in the sterile neutrino parameter space 
sin 2 20 ee — Am^j for the three current reactor experiments, namely Double Chooz, RENO 
[35] and Daya Bay. A similar study was performed in [28] for Double Chooz and Daya Bay, 
where a more approximate analysis in the sin 2 20 ee — Am^ plane was done assuming an 
overall systematic error and neglecting detector resolution. In our present work, we have 
used simulations with reduced values of errors, as quoted in the technical reports of the 
individual reactor experiments, where the cancellation of correlated reactor-related errors 
by using both near and far detectors has been taken into account. Also, we have considered 
realistic detector resolutions, which play an important part in the sensitivity analysis. The 
allowed regions in sin 2 20 ee — Am^ plane presented in this work are relevant in view of 
the fact that there is still significant tension in the existing data between appearance and 
disappearance experiments. Although an exclusion region is quoted in the present global 
best-fit scenario [11], we have tried to survey critically the region which would be accessible 
by this specific set of new reactor experiments, to see what it adds to existing information. 
We also compare our results with those of older reactor experiments like BUGEY [36], 
GOSGEN [37] and Krasnoyarsk [38]. 

The paper has been organized as follows. In Section 2, we briefly present the theory 
of neutrino oscillations with sterile neutrinos. In Section 3, we describe the details of the 
three reactor experiments we study, as listed in the technical reports of each experiment. 
We also give details of systematic errors and detector resolutions used in our analysis. In 
Section 4, we present our results, and conclude with a discussion of the results in Section 



2 Theory of Neutrino Oscillations with Sterile Neutrinos 

We know that neutrino mass and flavor eigenstates can be related by the relation 



where U is a unitary matrix. In matrix form, this mixing between flavor eigenstates v a {ct = 
e,fjL,r,s; where s stands for sterile neutrino) and mass eigenstates Uj(j = 1,2,3,4) in four 
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neutrino scenario can be represented as 



(u e \ /U el U e2 U e3 U e4 \ 

Vy, Ufja U^2 U^3 UpA 

U T U T l U t2 U t3 U T 4 

\V S ) \Usi U s2 U s3 U sA J 



This unitary matrix can be parametrized [39] as 



v 2 

V 3 
W 



(2.2) 



U = R 



34 ^34 



, 0)^24(^4, 0)R 2 3(9 2 3, 5 3 )Ri4(e 14 ,0)R 13 (6 13 , 5 2 )R 12 (9 12 , (2.3) 



where Rij(9ij,Sk) is the complex rotation matrix in the i-j plane and the elements are given 
by 



[R 
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or 
and 
and 
and 
otherwise 



q = j 
j 

i 



(2.4) 



Here, 9ij is the angle of rotation in i-j plane [17, 39]. If we assume #14, #24, #34 = then the 
above mixing matrix reduces to the standard Pontecorvo-Maki-Nakagawa-Sakata (PMNS) 
[40-42] form. Here we consider the (3+1) scenario, which is an extension of the three- 
neutrino scenario with the addition of one massive sterile neutrino. Solar and atmospheric 
neutrino analyses strongly discourage the (2+2) scenario due to the absence of sterile 
signals in the atmospheric parameters [43]. The (3+2) scenario is favoured by the tension 
between appearance and disappearance experiments in the (3+1) case, but disfavoured by 
cosmology. 

The three new reactor experiments we have considered have two detectors, one at the 
far site and the other at the near site. In these experiments, the electron antineutrino 
disappearance probability is measured through electron antineutrino events. The far de- 
tector is placed at sufficient distance from the reactor site (compared to the near one) so 
as to maximise the disappearance of the electron anti-neutrino. Such a two-detector set 
up has a great advantage over a single detector, as it can cancel or reduce the systematic 
uncertainties. 

In the presence of a sterile neutrino, the standard 3-flavor neutrino oscillation picture 
changes, and hence the survival probability must be modified due to the effect of (3+1) 
mixing. The baselines relevant to these experiments are of the order of a few hundred 
metres, and hence the oscillation would show signatures at the atmospheric scale as well as 
possible sterile-scale effects. The electron neutrino survival probability expression relevant 
for our analysis is [28] - 



1 - cos 4 9 U sin 2 29 13 sin 2 (^|^) - sin 2 29 14 sin 2 ( Am ' 4jL ' 



4E 



(2.5) 
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Reactor No 


Near detector(km) 


Far detector (km) 


1 


0.47 


1.12 


2 


0.35 


1.00 



Table 1. Core to detector distances in Double Chooz 



where L is the baseline and E is the neutrino energy. Here, P ee is seen to be a function 
of two mass-squared differences (Am 2 4 and Amf 3 ) and two mixing angles #13 and #14. 
When #14 — > 0, we recover the standard three- flavour probability expression in which the 
solar mass scale is neglected. 

3 Details of the experiments 

In this section, we present some technical details of the experiments for the sake of com- 
pleteness, collected from [25, 27, 35]. 

3.1 Double Chooz experiment 

The Double Chooz reactor experiment [25, 35] is designed to detect v e through the inverse 
j3 reaction- 

u e +p^e + + n (3.1) 

The anti-neutrino flux coming from the two nuclear cores of the Chooz power plant results 
from the (3 decay of the fission products of four main isotopes- 235 U, 239 Pu, 241 Pu and 
238 i7. The Double Chooz far detector is cylindrical in shape, having a radius of 115 cm 
and a height of 246 cm, and hence a volume of 10.3 m 3 . Both near and far detectors are 
identical inside the PMT support structure and the mass of each is about 10.16 tons. The 
reactor site contains two reactors, each producing a thermal power of 4.27 GW. The near 
detector is situated roughly at a distance of 410 m from the reactor site, while the far one 
is at a distance of 1067 m. The reactor core-to-detector distances are tabulated in table 1. 

In our analysis, we have considered an exposure time of three years with a 12% detector 
resolution. The x 2 function we have used is taken from the collaboration report. 

3.2 Daya Bay 

The Daya Bay neutrino experiment [27] works with reactor generated electron antineutrinos 
and uses a gadolinium (Gd) loaded liquid scintillator detector. It has two pair of reactors at 
Daya Bay and Ling Ao I, which generate 11.6 GW of power. One more reactor site at Ling 
Ao II is under construction. Daya Bay consists of three underground experimental halls, 
one far and two near, linked by horizontal tunnels. Eight identical cylindrical detectors are 
employed to measure the neutrino flux. The mass of each detector is about 20 tons. Four 
of these eight detectors are at the far zone while two detectors are kept in each near zone. 
The distance of the detectors from the reactor cores at the Daya Bay site is 363 m while 
this distance at the Ling Ao site is 481 m. The far detectors are at 1985 m and 1615 m 



- 4 - 



rVcdClUi 1NU 


iNtJclI LlcLcCLUi 1 Kill 1 


r dl U.C leC IUI ^KIll J 


1 


70 


1.52 


2 


0.48 


1.43 


3 


0.32 


1.39 


4 


0.32 


1.39 


5 


0.48 


1.43 


6 


0.70 


1.52 



Table 2. Core to detector distances in RENO 



Name of Exp 


Double Chooz 


Day a Bay 


RENO 


Location 


France 


China 


Korea 


No of Reactor cores 


2 


4 


6 


Total Power(GW^) 


8.7 


11.6 


16.4 


Baselines- near /Far (m) 


410/1067 


363(481)/1985(1615) 


292/1380 


Target mass (tons) 


10/10 


40 x 2/10 


16.1/16.1 


No of Detectors 


2 


2 


2 


Exposure(years) 


3 


3 


3 


Resolution(%) 


12 


12 


12 



Table 3. Details of the three reactor experiments 



respectively from the Daya Bay and Ling Ao reactor sites. We have used an exposure time 
of 3 years, with 12% resolution. 

3.3 RENO 

The RENO experiment [35] is designed to search for reactor antineutrino disappearance 
using two identical detectors. The set-up consists of a near detector roughly 292 m away 
from the reactor array center, while the far detector is about 1.4 km away from the reactor 
array center. This design of identical detectors at the two sites helps in reducing systematic 
errors. Each detector contains 16 tons of liquid scintillator which is doped with gadolinium. 
There are six reactor cores in RENO, for which the core to detector distances are listed in 
Table 2. 

The average total thermal power output of the six reactor cores is 16.4 GW^, with 
each reactor core generating about equal power. The energy of the antineutrinos in these 
experiments is in the range of 1 to 8 MeV. The resolution is 12%. 

We present the particulars of the three reactor experiments in Table 3. 

The systematic errors associated with each experimental set-up are listed in Table 4. 1 

1 Since we did not find specific values of the scaling and overall normalization errors of the Daya Bay 
experiment in the literature, we have assumed values similar to Double Chooz. 
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Name of Exp 


Double Chooz 


Day a Bay 


RENO 


Reactor correlated error(%) 


2.0 


2.0 


2.0 


Detector normalisation error (%) 


0.6 


0.5 


0.5 


Scaling or calibration error(%) 


0.5 


0.5 


0.1 


Overall normalization error(%) 


2.5 


2.5 


2.0 


Isotopic abundance error(%) 


2.0 


2.0 


2.0 



Table 4. Systematic errors associated with the three experiments 



4 Details of old experiments 

The details of the three old reactor experiments ( BUGEY [36], Gosgen [37] and Krasno- 
yarsk [38]) with which we compare our exclusion regions for sterile parameters are tabulated 
in Table 5. 



5 Results 

We have generated the 90% c.l. exclusion plots for sterile oscillations for all the three 
current experiments and compared the results with the old reactor experiments Bugey, 
Gosgen and Krasnoyarsk. The right side of each contour shows the no-oscillation region 
while the region left to the contours is the possible oscillation region. The results are found 
to be very sensitive to the values of systematic errors. In our calculations, we have used 
GLoBES [44-48] for simulating the experiments. The details of the statistical analysis are 
as follows: the total no of bins used in all three experiments are 62 and the width of the 
energy window is 1.8-8 MeV. The resolution used is 12 %. The uncertainty associated 
with the shape of neutrino energy spectrum for all the three experiments is 2 %. In our 
calculation we have used the x 2 function as defined in the GLoBES manual and used in 
standard GLoBES sensitivity analysis: 



x 2 = \^ V" - (1 + Qfl + ad)T d ,i) + ^ + <4l + <4i (5 1) 

i=l d=N,F ' R N F 



where Ojv,i, Of,% are the event rates for the i th bin in the near and far detectors, calculated 
for true values of oscillation parameters; 7^ j are the expected event rates for the i th bin 
in the near and far detector for the test parameter values; or, op, ajy are the uncertainties 
associated with the reactor flux and detector mass; and ctr,<jf,0n are the respective 
associated standard deviations. 

Fig. 1 shows the variation of the average bound on sin 2 29 ee from each of the three ex- 
periments Double Chooz, Daya Bay and RENO as a function of the bin-to-bin systematic 
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Baselines (m) 


15,40,95 


37.9,45.9,64.7 


57,57.6,231.4 


Target mass(~tons) 


1.67 


0.32 


0.4 


No of Detectors 


1 


1 


1 


Exposure(~ years) 


0.2 


0.39,0.56,0.98 


0.06 


Resolution(%) 


6 







Table 5. Details of old reactor experiments 



sin 2 26 pp Vs Bin-to-Bin Errors 



DAYABAY 
RENO 
Double Chooz 



0.0001 0.001 0.01 0.1 

Bin-to-bin error 



Figure 1 . Dependence of the average sin 2 29 ee bound on the correlated bin-to-bin error for Double 
Chooz, Daya Bay and RENO. 



© 
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error for a constant overall normalisation. We depict the behavior with respect to the cor- 
related bin-to-bin error since this is found to have the most significant effect on parameter 
sensitivities. The figure shows that the dependence of the sensitivity in the case of Daya 
Bay is slightly steeper, even for low values of the bin-to-bin systematics, than for the other 
two experiments. 

In our further analysis, we have used a reduced set of systematic errors as listed in 
Table 6, taking into account the partial cancellation of errors due to the presence of both 
near and far detectors, as documented in the experiment literature. 

We have included the changed flux normalizations given by the reactor flux anomaly 
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Reactor correlated error(%) 


0.06 


0.087 


0.5 


Detector normalisation error(%) 


0.06 


0.12 


0.5 


Scaling or calibration error(%) 


0.5 


0.5 


0.1 


Overall normalization error(%) 


0.5 


0.5 


0.5 


Isotopic abundance error(%) 


0.06 


0.087 


0.5 



Table 6. Set of reduced errors used in our calculation 
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Figure 2. Comparison of Double Chooz exclusion plot with other reactor experiments at 12% 
resolution using modified errors. 

[3] in our sensitivity analysis. We are performing a simulation and not a data analysis, and 
hence varying the flux normalizations has minimal effect on our results, since the relative 
normalizations simultaneously affect both the true and test spectra in the simulation and 
lead to a cancellation of their effect in the parameter sensitivity. For the same reason, 
leaving the flux normalization as a freely varying parameter does not have a major influence 
on our results, and therefore they may be taken to be indicative of the effect of spectral 
information only. 

In Fig. 2 and 3, we show the comparision of the exclusion plots for Double 
Chooz, RENO and daya Bay experiments with existing reactor experiments using modified 
errors. Our results for all the three experiments show better exclusion regions for sterile 
neutrino oscillation parameters than Gosgen and Krasnoyarsk in the range Am^O.01 to 
1 eV 2 . Our bounds are an improvement over Bugey in the range A?72-4 1 =0.01 to 0.05 eV 2 
but in the Am^ region above this, Bugey gives better bounds. From these curves, we see 
that sin 2 26 ee > 0.1 and sin 2 26 ee > 0.08 is excluded for sterile oscillation in the Am 2 1 =0.01 
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Exclusion contours in the sin 2 20 ee -Am 2 plane with reactor anti-v e events 
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Figure 3. Comparision of RENO and Daya Bay exclusion plots with other reactor experiments at 
12% resolution using modified errors. 

to leV 2 region for Double Chooz and RENO respectively. The Daya Bay exclusion bound 
in the region Am| 1 =0.1 to 1 eV 2 is found to be nearly sin 2 29 ee > 0.07. 

6 Conclusions 

From the above results, we can draw the following conclusions: 

• The sensitivity of reactor experiments like Double Chooz, Daya Bay and RENO to 
sterile oscillation parameters is significantly dependent on the systematic errors, de- 
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tector resolutions and uncertainties of each experiment. The dependence is especially 
strong on the correlated reactor-related errors and the normalization uncertainty. 

• Because of the multiple detectors and baselines in each of these experimental set-ups, 
it is possible to have partial cancellations of the experimental errors, especially of the 
correlated errors, which is beneficial in giving us better parameter sensitivities. 

• In an analysis with duly reduced values of errors, it is possible to obtain better 
bounds with these set-ups than those from many of the older reactor experiments, in 
spite of the fact that the relatively higher baselines in this case are less suited for a 
determination of oscillations in the Am 2 4 = 1 eV 2 range. For the latter reason, we 
find that we obtain better bounds in the low region 0.01 < Amy < 0.05 eV 2 . 

It may be noted that this region lies below the present best-fit range for sterile-scale 
oscillations. However, the tension between the sterile parameter bounds from different sets 
of experimental data like the appearance and disappearance experiments indicates that 
there is still significant uncertainty on the favored region. It is possible that the global fit 
regions may shift appreciably in the future. In view of this, these results are significant 
in showing that present reactor experiments like Double Chooz, Daya Bay and RENO 
may be able to give improved sterile oscillation bounds in specific regions of the parameter 
space. Such signatures may contribute in modifying the overall picture of sterile oscillation 
parameter sensitivity. 
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